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SUMMARY 


An investigation was carried out to study the elastic behavior 
after buckling of a rectangular plate reinforced with longitudinal 
stiffeners and subjected to compressive loads in a direction parallel 
to the stiffeners. Two possible buckling modes were investigated; 
namely, the plate may buckle around the stiffeners as nodal lines and 
the plate may buckle with the intermediate stiffeners as a unit. The 
edge stiffeners were assumed to have finite torsional rigidity, infinite 
bending rigidity peipendicular to the plate, and either infinite or zero 
bending rigidiiy in the plane of the plate. For the first buckling mode 
the calculated results were compared with experimental results and the 
agreement was found to be good. The method of analysis ixsed is a modi- 
fied variational procedure. Instead of solving three nonlinear partial 
differential eqviations simioltaneously, the transverse deflection w is 
assumed in the form of a fimction which satisfies the boundary condi- 
tions but contains \andetermined parameters. In terms of this assumed 
expression for w, the other displacement components u and v may be 
solved from the differential equations. Then the unknown parameters 
in w are determined from the condition that the potential energy of 
the system must be stationary. It was found that this modified varia- 
tional procedure will give much better results with the same amount of 
conputatlonal labor. 


INTRODUCTION 


Aircraft engineers are confronted with the difficult task of 
designing monocoque structiires to fulfill ail strength requirements and 
yet to remain within the limits prescribed by wei^t economy. To obtain 
an efficient monocoque structure, the engineer must have accxorate infor- 
mation regarding the stress distributioh in the structure. One of the 
most importEint component parts of the monocoque structure is the rein- 
forced plate under the action of compressive loads. 
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A reinforced plate considered in this report is a sheet reinforced 
with stiffeners in the direction of the applied loads and equally spaced 
dividing the sheet into panels (fig* 1(a)). Under the action of the 
applied con^jressive load, the plate may deflect around the stiffeners as 
nodal lines as in case 1 or the system may displace as a complete unit, 
as shown in case 2, the stiffeners howing with the plate (fig. 1(h)). 

In case 1, it is possible to consider the plate as one of the panels of 
the continuous plate attached to the stiffeners. In case 2, it is nec- 
essary to consider the stiffeners and the plate together. 

As mentioned previously, to obtain the apportionment of the applied 
loads to the plate and the stiffeners it is necessary to know the stress 
distribution in the elastic structure. Up to and including buckling, 
the stress distribution in the panels may be satisfactorily predicted 
by the linear- or small- deflection theory. After buckling, however, the 
lineax theory fails and the large- deflect ion theory of plates must be 
used. The large -deflection theory requires the solution of nonlinear 
differential equations which are generally intractable. In many cases, 
however, good approximate solutions can be obtained by the energy 
method or the variational method. 

In this report an investigation was carried out to determine the 
stress distribution after buckling in rectangular plates reinforced with 
longitudinal stiffeners and subjected to conipressive loads in a direc- 
tion paxallel to the stiffeners. Instead of following the energy or 
variational method as generally enployed, a modified procedure is used. 

It may be recalled that the fundamental differential equations to be 
solved consist of three simultaneous nonlinear equations in terms of 
displacements u, v, and w. However, if w is assumed to be of some 
fom, two of the differential, equations become linear in u and v 
and can then be solved. Therefore an alternative way of carrying out 
the energy method may be as follows. First assume w in terms of a 
series which satisfies the bo\mdary conditions but contains undetermined 
parameters. Substitution of w into two of the differential equations 
makes it possible to solve for u and v in terms of these undeter- 
mined parameters instead of by the usual way of assuming u and v in 
the form of series with additional parameters. The values of these 
parameters may now be ascertained so that the potential energy of the 
system is stationary. This condition implies that the third differential 
eqxiation is approximately satisfied. This alternative procedure has two 
advantages. First, it is relatively easy to represent accurately the 
deflection surface w by a series of a few terms because of the 
existing experimental observations. Second, without any additional 
parameters in u and v only the parameters in w need to be deter- 
mined in this case and much more accurate results can be obtained than 
could be obtained otherwise by the same amount of computational labor. 
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The behavior of a rectangular panel after buckling according to 
case 1 of the buckling mode has been treated by Dunn (reference l) . 

One may, however, take issue with Dunn’s work on two coimts. First, 
instead of satisfying the three equilibrium equations, Dunn's solution 
satisfies only two of the three equations and the third equation is 
replaced by an en^jirical relationship. Second, in determining the 
buckling load Dunn used one set of boundary conditions, but in studying 
the behavior after buckling another set of boundary conditions was used 
which was not at all compatible with the first set. 

In this report, analysis has been carried out for both case 1 and 
case 2 of the buckling modes. Deflection pattern, stress distribution, 
and effective width are calculated to con 5 >are with Dunn's results and 
the experimental results obtained by Ramberg, McPherson, and Levy 
(reference 2 ) . 

The work was carried out at the Daniel Guggenheim School of 
Aeronautics, College of Engineering, New York University, under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. The authors are indebted to Professor F. K. 
Teichmann, Chaiiman of the Department of Aeronautics, for his under- 
standing cooperation. 


SYMBOLS 


a 

Ast 

b 

C 

C^T 

d 

D 

®^®st 

G 

I 


length of rectangular plate 
stiffener area 

width of rectangular plate 

torsional stiffness of a sturdy stiffener when attached to 
rectangvilar plate 

torsion-bending constant of stiffener sectional area 

distance between stiffeners 

plate stiffness per unit width ^Et3/l2(l - 

Young's modulus for plate and stiffener, respectively 

modulus of rigidity 
moment of inertia 
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Ip polar moment of inertia of stiffener sectional area atout 

eixis of rotation 

J torsion factor of stiffener 

K tuckling constant 

m numter of half waves into which plate huckles in x-direction 

n,ng^ numher of hays and stiffeners, respectively 

p lateral pressure 

V strain energy 

u,v X- and y- components of displacements, respectively 

w lateral plate deflection 

Wg effective width 

x,y,z coordinate axes (fig. 2 ) 

P aspect factor (xA>) 

7 shearing strain 

S' = '1* 

axial strain in x- and y-direction, respectively 
constant of propoil tonality defining location of stiffener 
half-wave length (a/m) 

p' = dD/ 2 C 

V Poisson's ratio 
0 - 


X 

p = hD/C 


cr 


con5)ressive stress in stiffener 
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st 


T 


normal stress in x- and y-direction, respectively 
stiffener stress and strain, respectively 
shear stress 


cp = jrm/ 


a 


cr 


quantity at buckling, used as a subscript 


GOVERNING DIFFERENTIAL EQUATIONS AND BOUNDARY CONDITIONS 


In terms of the displacement components u, v, and w and con- 
sidering the case where w is large compared with u and v, the dif- 
ferential equations of equilibrium for a thin flat plate may be written 
as 


^x ^ 1 - V 1 - V ^x^xx ■*' ^y^xy) 


^x^yy " ° 


( 1 ) 


^xx + V (V V) l^(^x^xy + ^y^yy) + Vy " '^x^'xy = ° 


( 2 ) 


" D “ i ^X^)(^XX + VWyy) + (Vy + I ^y^ ) (Wyy VW^) + 

2(1 - V) (uy + v^ + Vy)^x^ = 0 (3) 


where 




is the lateral loading on the 


plate, and the subscripts denote partial differentiation. The median- 
fiber strains are 
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:« = +1 w/ 


6 ’ = V + W ^ 

y y 2 y 


7 *=U + V + WW 

'xy y X X y 


y 


j 


(^) 


and the median-flher stresses are 

E 


= 


E / • 

1 - 


+ VSy‘) 


« E 


a 


1 - V 
E 


2(v' ■" ''"x’) 


^xy 2(1 + v) ^xy J 


( 5 ) 


The extreme-fiher tending and shearing stresses are 


GTx = 


a = — 


Et / , \ 

^2j("XX + VWyyj 


Et / 


+ VW ^ 

yy xx; 


Et 


'"xy ■ 2(1 + v) 




( 6 ) 


Equations (l) to (3) niay he obtained hy substituting relations (. 4 ) 
and (5) into the three equilibrium eqiiations as given in reference 3 > 
page 305- 

To formulate the boundary conditions, it is necessary to survey 
briefly the general picture of the structure to be analyzed, the condi- 
tions which exist within it, and the loads that are imposed on it. The 
structure consists of a thin flat plate to which are continuously 
attached parallel and equidistant stiffeners. It is assumed that the 
length parallel to the stiffeners is leirge in comparison with the 
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stiffener spacing. The structure is loaded with a con 5 >ressive load 
parallel to the stiffeners in such a way that, the axial loads are 
identical in each stiffener. 


The stiffeners along the edges of the plate are assigned to have 
Infinite bending rigidity in the z-direction; that is, 

w = 0 at y = 0 and b ( 7 ) 


In the y-direction, the sizes of the stiffeners are usually such that 
the bending rigidity is much smaller than that in the z-direction. In 
the theoretical treatment, two limiting cases may.be taken. First, the 
stiffeners may be assumed to have zero bending rigidity. In such a 
case, the botmdary condition is that along the edges the median-fiber 
stress Cy is zero; that is, 

a =0 at y = 0 and b (8a) 

*7 

On the other hand, the stiffeners may be assumed to have infinite 
bending rigidity, and, in such a case, the boundary condition is 

V = 0 at y = 0 and b (8b) 


In the x-direction, it is assimied that the stiffeners are under a 
xmiform conipressive strain that is. 


u = -eg.j.x at y = 0 and b 


(9) 


The slopes of the deflection siirface along the ekges depend on the 
torsional rigidity of the stiffeners. Along the edge of the plate 
adjacent to the stiffener at y = 0 the internal bending moment may be 
written in the form; 


“y = '^('^yy ^xx) 

where D is the flexural stiffness of the plate. 


( 10 ) 


Assuming the plate cut along the stiffener, the moment given by 
equation (lO) may be considered as an externally applied moment acting 
on the stiffener. This bending moment, which acts at any point x, 
produces twisting in the stiffener. The moment per unit length of the 
edge is numerically equal to the rate of change of the twisting moments 
in the stiffener at the point x. Expressed analytically. 
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For a sturdy stiffener attached, to a plate in compression, it was 
shown by Lundquist and Stowell (reference 4) that the moment required 
to rotate the plate throu^ 1 radian is 


C 




( 12 ) 


where X is the half wave of the buckled pattern of the plate in the 
direction of the applied load, GJ is the torsional rigidity of the 
stiffener, cr is the uniformly distributed compressive stress in the 
stiffener. Ip is the_ polar moment of inertia of stiffener sectional 
area about axis of rotation, and (IgQ, is the torsion-bending constant 

of the stiffener sectional area about the axis of rotation at or near 
the edge of the plate. Thus 




where 9 is the angle of twist. Since 


(13) 


9 




8w 


the boixndary condition (U), aLfter substitution, becomes 

®(^yy ^ ^^xx) = -^xxy 

Similarly, along the edge y = b one obtains 

D(Wyy + VW^) = CWx^ (15) 

Consider the limiting cases of zero torsional stiffness and 
infinite torsional stiffness. In the case of zero torsional stiffness, 

C = 0 and equations (l4) and ( 15 ) reduce to My=0 at y=0 and b, 
which corresponds to the case usually referred to as simply supported 
edges. In the case of infinite torsional stiffness, C = » and equa- 
tions (l4) and ( 15 ) reduce to 0 = dw/dy =0 at y = 0 emd.b, which 
corresponds to the case us\ially referred to as clamped edges. 


METHOD OF SOLUTION 


Equations (l), (2), and ( 3 ) are three nonlinear partial differential 
equations to be solved simultaneously. The exact solution of these 
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equations is very diffiCTilt, if not impossible, witb the present 
boundary conditions. Appraximate solutions, however, can be obtained 
by means of a variational method or the finite-differences method. In 
the present work, a modified variational method is to be used. 

In the theory of elasticity, it is a well-known fact that the first 
variation of the potential energy of the elastic body with respect to 
the displacements u, v, and w leads to the differential, equations of 
equilibrium which are, in the present case, equations (l), ( 2 ), and (3) • 
The potential energy is defined as the difference between the strain 
energy and the virtual work which the surface stresses do over that 
portion of the boundary where these stresses are prescribed. For stable 
equilibrium, this stationary value is a minimum. The fact that the 
potential energy has a stationary value can be easily utilized in trans- 
forming the variational problem into an ordinary extreme-value problem 
by the Rayleigh-Eitz procedure as follows. Assume expressions for u, 

V, and w in terms of infinite series which satisfy the given boundary 
conditions but contain undetermined parameters. Substitute these 
expressions into the potential- energy expression and set the first 
derivatives of the resulting expression with respect to these parameters 
equal to zero. These parameters are then determined by solving the 
resultant simultaneous equations. With the values of the parameters so 
determined, the resulting expressions for u, v, and w satisfy both 
the bomdary conditions and the differential equations. If infinite 
terms of the series are retained and the series consist of con^jlete sets 
of functions, the solution may be proved to be exact. For the numerical 
determination of these parameters, however, only a finite number of 
terms can be retained. In such a case, the solution is then only 
approximate. 

In the case of a thin plate, the deflection surface form can be 
easily observed from ejqjerimental tests, and consequently the deflec- 
tion w can be represented fairly accurately by a series with only a 
few terms. It may also be observed that equations (l) and ( 2 ) reduce 
to linear differential equations if w is known. Therefore, for the 
same amoxmt of con5>utational labor, much more accurate results can be 
obtained by carrying out the variational procedure as follows. Assxime 
an expression for w in terms of a series which can represent well the 
deflection from observed tests. The series satisfies the boimdary con- 
ditions and contains undetermined parameters. Substitution of this 
assumed w into equations (l) and ( 2 ) results in two linear differential 
equations in u and v. Solving these equations simultaneously, "exact" 
solutions of u and v in terms of the parameters may be obtained. 

With the expression for w so assumed and expressions for u and v so 
obtained, these parameters can be determined by the usual Raylei^-Ritz 
procedure. Following this procedure, only a few parameters will be used 
and therefore there are only a few resulting simultaneous equations to 
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le solved. The amotmt of labor in the computation can thus he greatly- 
reduced. 

To satisfy the hounda^ conditions of \r, the deflection surface 
may he ass-umed "to he of the form: 

I = |^„(l - cos sj as (16) 

where Wq and are arbitrary deflection amplitudes -fco he deter- 

mined and t is the thickness of the plate. The form of equation (16) 
is selected so that w satisfies the conditions of simply supported 
edges when ¥q = 0 and the conditions of clamped edges when \I± = 0. 
The ratio therefore a measure of edge restraint. 

This ratio Wq/w^ is related to the elastic restraint of the 
sides through the houndary condition given hy equation (lU): 

D^Wyy +VW^^) = -CW^ (17) 

Since w = 0 at the side edges and w is assumed "to he of the form 
given hy equation (16), the foregoing equation may he written in the 
form 


DWyy = Ccp^Wy (18) 

at y = 0, where cp = inn/a. 

Substitution of equation (16) for w in equation (18) yields the 
following relationship 

where p = hD/c is the ratio of the flexural rigidity of the plate 
he-tween the edge stiffeners to the torsional rigidity of the stiffeners 
and P = \/h. 

The deflection- surface foim may now he given in terms of the 
maximum amplitude Wq, with the replacing of W^^ in equation (16) hy 
equation (l9)^ an<f is as follows; 
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^ = Wq (1 - cos 2 pcpy + A sin pcpy) sin cpjc 

X 


( 20 ) 


4 p 

where A = — (iP . 

It may he noted that at y = h the condition given hy equa- 
tion ( 15 ) is also satisfied. For plates with simply sipported edges p 
becomes infinite, Wq becomes zero, and then AWq becomes as 

defined before. 


Expressions of u and v 

The expressions for the displacement functions u and v may be 
obtained by solving the simultaneous differential equations (l) and (2) 
with the given w expression (20) . These equations, being linear in u 
and V, may be solved as follows. 

The complementary solution of equations (l) and (2) may be 
obtained by solving the following homogeneous equations; 


1 - V ^ ~ 


1 + V / . _ 

V + V + IvL + V \ = 0 

XX ''yy 1 - V yy) 


( 21 ) 




To fiilflU the boundary condition (7)^ "the appropriate form of the 
complementary solution can be shown to be as follows: 

u = (pt%^^ + Cgj cosh 2cpy + 2 C 2 <py slnh 2cpy + sinh 2qpy + 

2 Cij^qpy cosh 2qpj^ sin 2qpx - (22 

V = yt^o^ i ^ 2 C 2 - sinh 2cpy - 2C2Cpy cosh 2yy + 20^ - 

cosh 2qpy - 2Cij.cpy sinh 2yy j cos 2qpx + (2- 
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where cp^t^6g-t “■ ^2^ ^3-» ^6 inte- 

gration constants. Physically, the complementary solution gives the 
displacements due to tending. 

The particular Integrals of the differential eqviations (l) and (2) 
give the additional displacements u and v due to large deflection. 
The method for determining these particular integrals is given in 
appendix A where the following forms for u and v are obtained: 


u = q)t%Q^^EQ + E]^ sin Pcpy + E2 cos 23q)y + E^ sin 33<?y + 

E|^ cos U3<j)y) sin 2qix (24) 

V = cpt^Q^P cos Pqpy + F2 sin 2Pq^ + F^ cos 33<Py + 

F]^ sin 4p(py) cos 2<px + ^F^ cos + Fg sin 2pcpy + ¥j cos 33cpy + 

Fq sin 4pq)j^ (25) 


The expressions for the E's and F's are to he fo\md in appendix A. 

The general solutions for the displacement functions u and v 
are the sums of the complementary solutions and the partici0.ar integrals. 
The integration constants C2> and Cg may now he 

deteimined from the given houndary conditions as is carried out in 
appendix B. With the houndary conditions satisfied, the general solu- 
tions may now he written as: 


u 


cpt%Q^ +. C2 - 2C2cpyj(cosh 2q?y - sinh Pcpy) + Eq + E^ sin p(py + 
E2 cos 2Pq)y + E^ sin 3PTy + cos 4pcpy + E^ sin 539yl sin 2cpx - 


cp2t2rgtx 


(26) 
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V = 


cptV2. 


'Cn - 


1 - V 


1 1 + 


2C2 - 2C2<ryj{cosli 2cpy - sinh 2(py) + 


cos Pqpy + F^ sin 2P(py + F^ cos SPipy + sin 4pqfy^ cos 2cpx + 
P^C^ + CgCpPy + F^ cos p<py + Fg sin 2P(py + F^ cos 3Pcpy + 


F0 sin 4pcpy^ 


In equations (26) and (27), C^, C^, C^, and Cg are now known 

functions of the deflection parameters. The expressions for these 
coefficients are to he found in appendix B. 


(27) 


Evaluation of Potential Energy 


As mentioned previously, the potential energy is the difference 
between the strain energy and the virtual work. The total strain 
energy of the panel is the sum of the strain energies of the plate and 
stiffeners. The strain energy of the plate is composed of two parts, 
one due to the bending of the plate V-j^p and the other due to the 

extension or stretching of the plate Let denote the total 

plate strain energy, then 


where 


V = 
P 


hp 


+ V 


ep 


(28) 


Et- 


"" ■ 24(1 - v^) 



^^xx + ^yy)^ 


2(1 - v)(w^Wyy - w^2) dx dy 


(29) 


and 


^ep 


Et 


5(1 - v^) 



'f (V’f +|(l - v)(r^'f] dx dy 


(30) 
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The above integrations are to be taken over the total plate area. For 
rectanguleir plates with edges held in the original plane it can be 
proved that 





w- 


xy 


■) dx dy 


= 0 


(31) 


With the assigned expression for w in equation ( 20 ), and the 
expressions for u and v given by equations (26) and (27)^ integra- 
tion of equations (29) and (30) gives 


^ EyVab 2 
2 1 o 

1 - V 


(32) 


and 


V = + «3"stWo^ + (33) 

The evaluation of these expressions is carried out in detail in 
appendix C. The' expressions for Mg, M^, and are given by 

equations (C 3 ) ^ (C 6 ), and (Cj) in appendix C. 

The strain energy of a stiffener Vg is conrposed of three parts, 
namely, the extensional strain energy Vgg, the bending strain energy V-j^g, 
and the twisting strain energy The total strain energy is 

therefore 


^s ^es ^bs ^ts 


( 34 ) 


and the expressions for V^g, Vj^g, and V^g are: 


^es 2 ^st^st ^ ]P y=Tib ^ 


( 35 ) 


1 pQ- 2 

^bs " 2 ®st^st ("xx) y=qb ^ 


(36) 
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(37) 


where Tih gives the location of the stiffeners t] = ^/n, n is the 
number of hays, and j is an integer which is equal to or less than n. 
The evaluation of equations (35)> (36), and (37) is also carried out in 
appendix C. After some simplification, it is found that 


Ves = I 



+ Nr 


2W, 


) ^2^st 


est^ 


(38) 


'hs 




(39) 


V 


ts 




m 


The expressions for N^, N^, N^, and N|^ 

Introducing the following parameters 

^ ®st^st 

^ =-EdT 

g, ^ ^st^st(^ 
Edt 


are given in appendix C. 



P 


®st^st 

dD 




2C 
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where d = h/n is the distance between stiffeners, the strain energy 
of A stiffener then becomes 

^ * I; ?Bt^] 

(41) 

The virtueil work W is the product of the external load and the 
displacement in the direction of the load over that portion of the 
boundary where the boundary load is prescribed. In the present case, 

W is 


^A_j_ + Ebt^ 


= ^st ^ 


(^st’^st" st 


(^2) 


where n^^ is the number of stiffeners. 


The potential energy U of the panel is the simi of the strain 
energy of the plate and the strain energy of the stiffeners minus the 
virtual work 


Ij. 5 

y _ Eqp t ab 
1 - 




/ P^3 

TT 2 

^1 + -ii55— + m^i 

Wo + 


Ust^ ’ 


/ 2\^^st -^st 

^ “st + ij 




m 


In writing equation (43) it is assumed that all the stiffeners have the 


same properties 


, and ^ i 


indicates the sum of Nj 


T]b = d, 2d, 3d, and so forth. 


with 
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Determination of Unknown Parameters 

In the previous pareigraphs, the solution of the differential 
equations (l), (2)^ and (3) has been carried out in the following 
manner. Assume an expression for w which satisfies the boundary 
conditions but contains undetermined parameters. Solve equations (l) 
eind (2) in terms of the assumed w expression. If this assmed w 
also satisfies equation (3)^ then the solution is exact. In the present 
einalysis an atten^jt will be made to satisfy equation (3) approximately 
by using the condition that at equilibrium the potential energy has a 
stationary value. In reference 1 a similar problem was treated by Dunn; 
the difference between Dunn's method and the method used in this report 
is that Dunn determined the parameters in w by -empirical relations 
instead of by equation (3) • 

The parameters to be detennined in w are Wq and p. Inspection 
of equation (43) shows -that the po-tential energy U is of the form 

u = - (Rp^st - ®3K^ + I^4"st^ 

where and R^ are functions of 3 and R^^ is a cons-tant. 

It is observed from experiments that once a plate buckles, the 
number of buckled waves remains -the same when the buckling strain is 
exceeded. The factor 3, which defines the number of buckles in the 
direction of loading, is therefore to be taken as a constant after 
buckling and equal -to the value at buckling. 

The parameter Wq may now be determined from the condition that 
the potential energy U has a s-tationary -value with respect to Wq* 

The condition = 0 leads -to 

>^o|^lWo^ - {^2^ St - R3)] =0 (^ 5 ) 

Equation (45) gives three roots for Wq, namely 

Wq = 0 (46) 


and 


Wo = ±1 


l^^st " ^3 
2R, 


(4?) 
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The root Wq = 0 represents the unhuckled state. The positive and 
negative signs in equation (47) indicate that the plate may buckle 
indifferently in or out. Since these two possibilities eire Irrelevant 
to the present analysis, the positive sign will he taken in the sub- 
sequent discussions. Equation (4-7) also shows that when 


^st 



there is no real root and therefore ¥q = 0 is the only possible 
solution for equation (45). Physically, it indicates that when e gt 
is below a ceidiain value, the panel will not buckle and the unbuckled 
state is the only possible equilibrium form. The limiting condition 
occurs when 



m 


which indicates the buckling of the panel. 

Substituting equation (48) into equation (47), one obtains 



(49) 


The buckling stress of a plate is usually written in the following 

fonn: 


a 


cr 


12(1 - v^)b^ 


(50) 


where K is a nondimens ional coefficient that depends on conditions of 
edge restraint and shape of the plate. The buckling strain e is of 
course related to Ucr according to Hooke’s law. Rewrite equation (5) 
as follows: 


+ V0 


•) 


e("^x y 




(51) 
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When the houndary conditions are that Oy' = 0 along y = 0 and y = h, 
at huckllng stress a ’ = 0 everywhere. Therefore, one obtains 

J 

^cr “ ^ ct /^ 


£ind 



Hence 


p2 Eg 


(52) 


On the other hand, when the houndary conditions are that v = 0 
along y = 0 and y = h, at buckling stress it can he seen from equa- 
tion (b 4) that = 0 everywhere. From equation (5l) cTy’ = 
and 6^' = ^1 - v^jcr^yE. Therefore, one obtains 


^cr 




and 

K . jg L- -V ! l h ' (53) 

K 2 

Equations (52) and (53) give K as a function of p = a/mh. For a 
given value of m, K may he plotted against a/b for a given value 
of p which specifies the edge-support conditions. In a given problem, 
when a/h and p are known, the value of K may he obtained by 
choosing the mlnimimi value from the curves. The corresponding value 
of m is the number of half waves which actually occur. In practical 
cases, the dimensions of the panels are usually such that the curves K 
against a/b become strai^t lines, and ■ the value of K approaches the 
minimum value. For these cases, the values of p may be' taken to be 
those which correspond to the minimum values of K. With P determined, 
Wq may be found from equation (49) at any 6gf The median-fiber 
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strains and stresses may "be computed from equations (4) and (5) and 
the extreme-fiber tending stresses, from equations (6). 


RESULTS AND DISCUSSION 

Case 1: Rectangular Plate with Edge Stiffeners 


As mentioned in the introduction, when a stiffened rectangular 
panel is under the action of con^pressive load two huckling modes may 
occur. One case is that the plate may deflect around the stiffeners as 
nodal lines. In such a case, it is possible to consider only one of 
the panels between two stiffeners and to regard it as a rectangular 
plate elastically st: 5 )ported along the two unloaded edges by stiffeners 
having infinite bending rigidity in the z-direction but finite torsional 
rigidity.^ 


Neglecting the extensional strain energy of the stiffeners and 
letting b in this case be the width of the panel, the potential 
energy of the panel becomes 


U 


r 


+ 

0 

Jl- 

( _ 
^Ms€st 

N- 

)/®st ^St 


E bt 

“st + 


•st 


M4 - 


(54) 


The condition 5 u/3Nq = 0 becomes 
At buckling strain (^stj^ 


(55) 


/cr 


M3(est)cr + % + (n4/24p) = 0 


(56) 


In terms of K as defined by equations (52) and (53 ) > if tte edge 
stiffeners are assumed to have negligible bending stiffness in the 
y-dlrection, that is, Uy’ = 0 at y = 0 and y = b. 
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0.1{053P^) + Iip^(l.o8o8 + 2.l6i5p2 + 
^P^)]yp^(o.iK)53pV^ + 1.08073^ + 0.75) 

(57) 

and if the edge stiffeners are assumed to have infinite bending stiff- 
ness in the y-direction, that is, v = 0 at y = 0 and y = b, 

^ pMi + + 

A /p"|(3 . tap") . ^ p2(l . v?2)m . A . vp")^"] 

(58) 

For given values of m- and P, K for these cases can be calculated 
from equations (57) or (58) and is plotted against P for various 
values of |j. in figure 3* For a large a/b ratio, it is pointed out 
in a previous paragraph that the minimum values of K or Kmin iiiay be 
used. The values of Kmin and the corresponding values of P are 
plotted in figures 4 and respectively. The values of P corre- 
sponding to may be calculated from the condition dK/dP = 0, and 

may then be con^juted. The condition dFC/dp = 0, however, results 

in an equation with twelfth power in P eind is difficult to solve. 

But with the approximate values of p obtained from figure 3^ one may 
use Newton's method to ceilculate more-accurate values. 


K = 1^^3^(0.4053 + o.8io6p^ + 
2.70193^) + (0.75 + 23^ + 


Now for a given p the values of p corresponding to may 

be determined from figure 5^ and Wq may be calculated from equation (49) 

as a function of In figure 6 the values of ^oXj(^ stj^ cr) " ^ 

are plotted against p. ' ^ 


With Wq and p calculated, it is thus possible to conqpute the 
stress distribution and the effective width of the plate at a given 
value of . 
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Substituting in equation (5) the expressions of e^' and gy' 

given hy equations (B3) and (b 4) in appendix B, the median-fiber 
stress cTjj' may he written in the following form: 


a ' 

X 


^ - v) (c^ - 2 qC^y + C2)(cosh 2cpy - 

sinh 2qy) + + vHq^ + sin gcpy + (g^ + cos 2Pqpy + 

(G3 + VH3) sin cos 2q»c + (vp^Cg + Gj^. + + 

(g^ + vH^) sin ^(py + ^Gg + VH5) cos 2Pq)y + ^Gy + VH^) sin 3Pcpy + 


(gs + VH3 ) cos ^3<py 


> - 


J 

1 - V" 


cp^t^r, 


st 


(59) 


The distributions of a-^' over the loaded edges of the plate with 
different boundary conditions are plotted in figure 7 as compared with 
those calculated by Dunn. 


The total load carried by the plate for a given stiffener stress 
is 


P 



<3y 


The effective width of a plate in coii5)ression is by definition 


P = 


-2wgta 


st 


where the negative sign indicates that is a compressive stress. 

The ratio of the effective width to the plate width may now he written 
as 


w 


_e 

h 





(60) 


It is shown in reference 1 that the distribution of stress at the 
nodal lines gives the actual load carried by the plate and any variation 
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of the stress dlstrihution between the nodal lines is balanced by the 
shearing stresses at the stiffener. Therefore, the effective width of 
the plate may be found f2rom equation (60) when cpx is tahen to be 
equal to Jt, 2n, and so forth; carrying out the integration and taking 
qac = 7t, equation (60) becomes 


w 


e 

b 


E 


(1 - v^) 


2\ E 


st 


“ 6 [29b (*^1 1 + V ^2) + G4 + VH4 + 


(^1 ^1 *^5 ■*" 3it (^3 ^^^3 ^7 ^ 


The ratio of total effective width to initial width or 2we/b is 
plotted in figure 8 for the two limiting conditions p = 0 and p = 00. 

To con5)are these results with test results given in reference 2 
where the panels were stiffened by z- shaped stiffeners, a panel with 
the following dimensions is chosen as a typical specimen: 

p = 0.198 


a = 19 inches 
b = 4.0 inches 
t = 0.025 inch 
P = 6800 pounds 
Ag^ =0.39 square inch 


With these dimensions and the boundeiry conditions ffy' = 0 one 
obtains P = 0.688 from figizre 5* The calculated length of buckle 
developed is X = 0.688 x 4 = 2.752 inches. The corresponding number 
of buckles in the direction of the length of panel is m = 19/2.752 = 6. 
or 7 buckles since m must be an integer. Observed data showed the 
existence of seven buckles with a corresponding measured length of 
buckle of X = 2.7O inches. The critical buckling stress may be calcu- 
lated utilizing the coefficient K obtained from figure 4, and is 
found to be = 2570 psi, to compare with the observed value of 

= 2500 psi. 

The wave form of the plate between the stiffeners under the com- 
pression load of 6800 pounds may be obtained in the following manner. 
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For n = 0.198 and 3 = 0 . 688 , the amplitude factor 

V('stAcr) - 1 

is 0.538. The stiffener strain corresponding to the test load of 
6800 pounds is = 11 X 10 "^. Using equation ( 20 ) , the results of 

the calculation of the wave form throu^ the crest of a huckle x = \f 2 
are shown in figure 9 - The wave forms calculated hy the present method 
show good agreement with those observed from the test. 

With the hoimdary conditions v = 0 , the value of 3 from fig- 
ure 5 is 0.835. The calculated length of huckle developed is 
X = 0.835 X 4 = 3.3^ inches. The corresponding number of buckles in 
the direction of the length of panel is m = 19 / 3 - 3 ^ = 5*70 or 
6 buckles since m must be an integer. This is compared with the test 
observance of seven buckles having a measured length of \ = 2.70 inches. 
The critical buckling stress, utilizing figure 4 for the value of K, 
is calculated to be 2170 psi, as compared with the observed data 
of = 2500 psi. 

The wave form of the plate between the stiffeners tinder the com- 
pression load of 6800 pounds is calculated in a similar manner as for 
the boundary condition cr^' = 0 utilizing figure 6 for the amplitude 
W ^ 

factor . and is shown graphically in figure 9* 

V(^st/^cr) “ ^ 


The experimental values of 2 UeA> from these tests are also 
plotted in figure 8 for conparison. It is seen from the figure that 
with boundary conditions Uy' =0, the calculated values give very good 
comparison. It thus seems that, for the z- shaped stringers used, the 
boundary condition can be well approximated by assuming dy' =0 along 
the stiffeners. 

It may be pointed out here that in reference 5 Dunn implicitly 
used the condition dy’ = 0 in the determination of the buckling stress 
and then he took the results obtained in this case together with the 
boundary condition v = 0 in the study of the behavior after buckling. 
It is easily seen that these two boundary conditions are not compatible 
at all with each other. 


Case 2 : Rectangular Panel with Intermediate Stiffeners 

When the longitudinal stiffeners have small bending and torsional 
rigidities, the rectanguleir panel will buckle as a unit with inter- 
mediate stiffeners bowing out with the plate. The edge stiffeners are 
ass\mied to have infinite bending rigidity in the z-direction. Pour cases 
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will be considered, namely, the edge stiffeners having either zero or 
Infinite bending rigidity in the y-direction, together with either zero 
or Infinite torsional rigidity. 

Consider first the case where the edge stiffeners have zero 
bending rigidity in the y-dlrection and zeir> torsional rlgldityj |i = 00 
in this case and AWq becomes W-j^. In terms of W, the potential 
energy U of the panel in equation ( 43 ) becomes now 
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The condition Bu/SWi = 0 leads to 



The buckling strain is therefore 



From equation (64), one obtains 



With these equations, one may proceed with the cong)utation as in 
case 1 when the number of intermediate stiffeners is given. 

Consider next the case where the edge stiffeners have zero bending 
rigidity in the y-direction and infinite torsional rigidity. The 
potential energy U of the panel is given by equation (43) . With 
p = 0, one obtains 
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= [2 + (1 + I96 

- ^ 

2 = [35(1 - V^)/256 + |(l - v^)^^[2 + (1 + P^) ^^32^- 

j^(l - V^)p^ (l + p^) ^[7 + (17 + 8v)3^ + 13P^ + 3P^y6W 


M = -3(1 - v^)/8 


M,. 


= (1 - V^)/2 


= P^-* 2 (1 + 2v) + vP^ - cpiilD^l + P^)(l + v)(cosh 2q)i]'b - 


sinh 2<pr\h) + ^ + P ^ i p^ ~ v) 
2 / 2^2 
2[1 + r) 


cos 2tT(\ 


]A(h- 


J/ 




3 1 1 

^2 ~ Q ~ 2 ^ 


= sin^KTi/4 
K|^ = 4p^Bin^2iai 




y ( 68 ) 




Then K, and Wq are given by ihe following expressions: 




Ml + j 

t[- 

!I3 + (n 4 /p)J i(8^ 

- 

M3 + 

(z®V”) 


K = 


(1 - v^) 

4p^ 


^ 48nMi + ^ [^N3 + 


(69) 


-nM^ + 


I 


(70) 


6Nc 
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wdmi + yl r«3 + K/**)] frst 

1 ^ — r o\ 1 **=^ 

96nM2 + 246* 2 (%^ ■*■ ^2 y 



(71) 


The third case to he considered is that of the edge stiffeners 
having infinite bending rigidity in the y-direction and zero torsional 
rigidity. The potential energy- U of the panel is given by equa- 
tion (62) with 


M2’ = 



- v^) (1 + 


+ 


1 - V 
3 - V 



256 


M^’ = -(1 + vP ^)/8 


r (72) 


Nn 


I 



^ ^ 2 <pnb)(cosh 2 q)qb 

3 - V / 


sinh 2 cpr^b) / 8 




J 

Mjp N2 *j and N]j^’ being the same as given by equation (63). 

One can calculate and Wj from equations (65) and (67) "by the 

use of new M2', M^', and N^’ as given by equations ( 72 ). From the 

definition of K in this case 


K = 


n2 ^nMi’ + [pFcj’ + (Nh’/ki’) 


4 p^ 




5’Wr 


( 73 ) 


The fourth case is that the edge stiffeners have infinite bending 
rigidity in the y-direction and infinite torsional rigidity. The 
potential energy U of the panel is again given by equation ( 4 - 3 ) with 
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{35^iaii9+^+<i|3- 


“2 




J 


U32 + 


1 


(1 - v)^(l + 2V + vP 
l6n(3 - v)(l + p2)^ 

= -(3 + i^v32)/8 


y (7^) 


N, = P ^ ^ ^ .[(2q)Tib - 3)(costL acpTib - sinh 2qnib)l + 

^ 2(1 + 32)2(3 - V) ^ 


" 2(1 . 3")" 


y 


^If M4, N2, and N4 being the same as given by equations ( 68 ) 

Eqmtions (69) and (7I) give (^st)cj. while K is 


K = 


^2 ^ |^N3 + (%/l^’)] 




-nM^ + 


Hb-»£ 


(75) 


Numerical calculations have been carried out in the case of one 
stiffdner in the middle of the plate. This case is important because 
in the case of multi stiffeners approximate analysis can be made by the 
substitution of the single -stiffener method. In this case, n = 2. 

The bending rigidity ratio p is plotted against Km-in in figure 10 

with 5 or 5’ equal to 0 and 0.4; the cor responding 3 is plotted 
in figure 11. In figure 12 are plotted Wo//(e st/ ^ cr) - 1 and 

WjJ \/(^8t/^cr) “ 1* Finally, the distribution of a^' along the loaded 
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edge of the plate is plotted in figure 13 at estAcr boundary 

condition Uy’ = 0 and in figure l4 at ^st/^cx = 10 for boundary 
condition v = 0- 


Daniel Guggenheim School of Aeronautics 
College of Engineering 

New York University 

New York, N. Y., December 29, 19^9 
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APPENDIX A 

DETERMINATION OF PARTICULAR INTEGRAIS 
FOR EQUATIONS (l) AND (2) 


The nonhomogeneous differential equations are: 


+ ^xy) = 


(Al) 


"'^xx 1 _ V ^ ^yy) 


where 


+ ^y^xy) - (Vyy " ^y^xy) 

2 

^ ^ (^x^xy ■*■ ^y^yy) “ ('^xx'”^y ” ^x'^^xy) 


(A2) 


^ (A3) 


Substituting the expression for w, equation (20)^ into equations (A3)> 
^(w) and |(w) are obtained as follows: 

« 

- ((3 - ^vp2) + (1 - vp2)^^ + 

+ (3 - Uv)p^A sin pqy - ijX + (9 - v)3^A sin 33<py - 
|[8 + 8(1 - V)p2 + (2 + 2 p 2 )a^ cos 2pcpy + 
i(l + ^p2) cos 43<py^ sin 2cpQC 
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iCw) 


cp3t2wo2p 

1 - V 

^(3 + 6p^)a cos 3P<py ~ ^ ^ ^ sin 2p(py + 

(1 + sin 4pq)^ cos 2q)x - •^(3v + 2 P^)a cos pcpy + 

^(3v - 6P^)a cos 3pcpy - ^2v + - sin 2P(py + 


-^(3 - 2 p 2 )a cos pcpy + 


(v - sin 4pcpy|' 

Substituting the expressions just obtained into equations (a 1) 
and (A2) for ^(w) and i(v) , it is possible now to solve the differ- 
ential equations and to obtain the particular Integrals of the displace- 
ment functions u and v in terms of the undetermined deflection 
parameter Wq. The method of solution used here is the so-called method 
of undetermined coefficients. Examining the expressions for ^(w) 
and |(w ) , the particular integrals for u and v may be assumed to 
be of the following form: 


u = cpt%^Q^(EQ + El sin Pcpy + E 2 cos 2pcpy + E^ sin 3Pcpy + 

Elj. cos 4pcpy) sin 2qpx (a4) 


= cpt%o^P |(^1 cos Pqjy + F£ sin 2Pcpiy + cos 3P<py 


Fjj. sin 4pqpy) cos 2q)x + (f^ cos Pq)y + Fg sin 2p(py + 


F^ cos 3P<py + F 0 sin 


(A5) 
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Substituting equations (A^) and (A5) i^'to equations (Al) 
and equating the coefficients of like terms, the E and F 
may be determined to be as follows: 


Eq =-^[(3 - ^Vp 2 ) + (1 - vP^)a^ 




+ 2(6 - Uv)3^_+ 


E2 = 


Eq = 


1 Jitli + (2 - v)e^ + 2^5 + (1 + 

+ r2)2 L 


16(1 + 


4(4 + 


[s + 2(18 - v)P^ + 


^^ = -T6 


F-, = 


i [12 - (16 + llv)P^ - 2p^A 

„o\ 2 L- -• 


l<-(4 + p2) 


^2 - , ^ 

l6(l+p2)2L 

r, 3 + (16 _ v)p2 + isp^A 


F4=-| 


F5 = 7^(3 v + 2 P^)a 
F^ = — i-[4v + (v - P^)a^ 

l6p2 *- 

F? = V - 2 p 2 )A 

‘ 12^2 

^8 = - 

32p2 


and (A2) 
coefficients 
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APPENDIX B 

DETERMINATION OF INTEGRATION CONSTANTS 
OF EQUATIONS (l) AND ( 2 ) 


The general solirfcions for the displacements u and v are the 
sums of the complementary solutions and the particular integrals: 

u = cpt^W^^ljCi + Cgj cosh 2qy + 2C^py sinh 2qy + + C|^j sinh 2cpy + 

2Ci(^qy cosh 2<py + Eq + E^ sin Scpy.+ Eg cos 2Pcfy + E^ sin 33qy + 

/ 

Ei^ cos sin 2cpx - (B1) 


V = <pt%Q 




' 1 - V 
VI + V 


2 C. 



sinh 2qy - cosh 2cpy + 


“ 2C4,c(y sinh 2qy + cos Pqy + 


sin 


2P<py + F^ cos 3Pcpy + Fi^ sin ^tPqy) cos 2cpx + P^C^ + CgcpPy 


F^ cos Pqpy + Fg sin 2Pcpy + Fy cos 3P<py + F3 sin 4pcpyjj- (B2) 

Consider first the case where the edge stiffeners have negligible 
tending stiffness in the y-direction; that is, 

Utt' = + veY*) = 0 along the edges y = 0 and y = b. 

X H „2 V y ^ ' 
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Ttie strains in the median surface of the plate are given hy equa- 
tions (i^). Substitution of the expressions for u, v, and v given 
hy equations (Bl) , (B2), and (20) into equations (4) gives the following 
strain expressions: 


e^' = 2 (c3_ + C2 + 2 cpCi^y) cosh 2 qjy + 2(03 + + 


2 <FC 2 y) sinh 2cpy + Gq + Gi sin Pqy + G 2 cos 2^cpy + 


G^ sin 3P(fy cos 2cpx + G],. + G^ sin P.qy + G5 cos 2Pcpy + 


Gy sin 3Pcpy + G3 cos 4Bqfy|- - 


€y' = (p%\‘ 


^ (ri ~' v" *^2 " *^1 " cosh 2qy + 


(B3) 


- <=3 


2cpC2yj 


sinh 2cpy + Hq + sin pcpy + 


Hg cos 2Bcpy + sin 3PfO^ cos 2cpx + P%5 ^ 

sin Pcpy + cos 2Pcpy + Hy sin 3Pcpy + H3 cos 4p(pyX (b 4) 


7xy' = 2CRJT 




cosh 2cpy + cos, Pqpy + Jg sin 2Pcpy +. 
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■where 


Gq = I + A^) 




Go = - 


?^(v - 


>(l + ^®) 


“3 = - 


4-(l|. + 9p2)2 


P^a(^v - 9P^) 


\ = |(3 + 


G - 3A 
^5 - T 


% = - + ^ 2 ) 


r - A 

^7 - - u: 


Gs = - 
^ 8 


Hq = - i P^^ + A^) 


8 


Hi ^ P^a(j<- - vP^) 


11 


k{k + p 2 )‘ 
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H2 = 




2(1 + 3^)2 


Ho = 


4(4 + 93^) ‘ 


% = I + a2) 


He = - M: 


=6 = 5 '-('* + 


E^=ivA 


%= -gv 


1 = ^ „\2 + '') 


{k + 3^ 


■^2 = 7 — 

U + p^) 


Jo = - 7 — TTT 33^A(1 + V) 

^ U + 93^)" 
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The hauiidary condition (8a) leads to the following three equations: 


(Hq + vGq) + (h 2 + VGg) + 2(1 - v)(^^ C2 - Cl) = 0 

+ K (% - + (% - ''^st = 

^Hq + vGq^ + ^H2 + VG2) + 2(1 - v) ^1 C2 - Cl - 2qCi^h) 

" ^3 “ 2cpb = 0 


0 


Kb6) 


cosh 2cph + 




These equations result from the stipulation that the boundary con- 
dition must be satisfied for all values of x; thus it is necessary that 
the coefficients of each trigonometric term in x vanish at the 
boimdaries. 

The boundary condition ( 9 ) that u = -6g^x at y = 0 and y = b 
gives the following relations: 


Cl + C 2 + Eq + Eg + E]^ — 0 

^Ci + C2 + 2C|j^cpb) cosh 2(pb + + 2CgCpb) slnh 2cpb + 






(B7) 


Eq + E2 + Ei^. = 0 


Equations (b 6 ) and (B 7 ) consist of five equations for the detenninatlon 
of the six integration constajits . One more equation therefore is needed 
in order to determine these constants uniquely. This equation may be 
obtained by specifying that v is a constant at a certain point. How- 
ever, since the actual value of v is immaterial to the problem, 
equations (b 6 ) and (B 7 ) -will be used to determine Ci, 09^ C^^ Ci^., 

and Cg, -with C5 undetermined. Solving these equations, one obtains: 
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Cl - + VGq) + (H2 + VGg) - (Gq + G2 + 



+ Gg + Gg 

^)] 




C2 

= - 4 ^[( 

vGq] 

I + (h2 + VG2) 

+ (1 - 

v) ^G q + G2 + 


Gi^ + Gg + Gg 





C3 = 

= -.aVv,[(: 

Hq + vGq) 

+ (H2 + VG2) 

A - 

1 + 

v)2 / 

V ro + ^2 


G4 + Gg + Gg' 

J cosh 2<pb - 1 
J sinh 29b 



Ci,= 


+ VGq) 4 


(1 - v) 

(Gq + G2 + 


Gij + Og + Gg) 

1 cosh 29b - 1 
J sinh 29b 


/ 

C6 = 

’ sF'’ '1 - (' 

+ vgJ 

- 




Since 2cpb = 2it/p and p = a/mb, for usual panel dimension e^^^ 
is large compared with e-S^ti and 1. It is therefore sufficiently 
accurate to take 


sinh 29b cosh 2cpb 

r; — 

2 ® 


(B9) 
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¥itli these simplifications, it is obvious that 


-C- 




-Cij. 


(BIO) 


Wow consider the case where the edge stiffeners have infinite 
bending rigidity in the y-direction. The boundary condition that v = 0 
at y = 0 and y = b then yields the following relations: 




^2C^-C3 + (Fi+F3)p=0 


= 0 


(rrr ^^2 - - SC^cpb) sinh 2cpb + 2Ci, - C3 


> (Bll) 


2 C 2 cpb^ cosh 2cpb - + F^^ = 0 


+ rtiJg - F^ - Fy — 0 


Solving simultaneously equations (Bll) and equations (B7) which 
result from the boundary condition u = -Cg^x at y = 0 and y = b, 

one obtains 
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" [(®0 + (r + V ^ + 3 (^1 + F3 j (1 

sinli 2q)b - 2cpbj 


COBb 


Cg - 


(Eq + Eg + Ej^) sinb 2(pb - + F3) (l + 

cosh 29b) j sinh 2cpb - 2cpbj 

^3 " [i + V ^(®0 + Eg + Ei^)(coah 2cpb - l) + p^Fj + F 3 ^(siiih 2q)b - 
2-pb] j (f Blhh 2qpb - 2q)b^ 


> (B12) 


- I^Eq + Eg + Eij^j^cosh 29b - 1) - P^F^ + 
F3) sinh 29b j ^ sinh 29b - 29b) 


= -(F5 + Fy) 

''6 = 1(^5 ^7) 

With the simplification given by equation (B 9 ), the values of Cj, 
^ 2 ) C4 become 


-C3 


Cg % -Cjj^ % 


Eq + Eg + Ej^^^ + 


1 + V 
3 - V 


- _(®0 + ®2 + ^ 4 ) - p ( Ei + F3I 


(Fi + F3)] 


(B13) 


Since C]^ ~ -C3 and C2 ~ -Cj^. in both cases, the general solutions for 
u and V may be written in the siigjler form as given by equations (26) 

and (27). 
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APEEiroiX C 


EVALUATION OF STRAIN ENERGY 


The bending strain energy of the plate is given by equation ( 29 ) 
and is 




Since 

Wxx = (P^Wq (-1 + cos 2 Scpy - A sin Pcpy) sin cpx 
Wyy = q)^WQ( 4 p^cos 2 ^q)y - AP^ sin 3 q^y) sin cpx 

by substituting in expression (Cl) and integrating over the area of the 
plate, the bending strain energy is found to be 


„ i, O.U ,, ,, £ 

V-bp = ^ MiWo 

1 - V 


vhere 


Mp = ^ 6 + 3(1 + + 3(1 + 3 ^)^A^ + ^(1 + (C3) 

The extenslonal strain energy of the plate is given by equation (30) 
and is 

“ 0 ( 1 "! I f {‘y’f * 

" (c4) 

The strains in the median surface of the plate are given by equa- 
tions (B3), (B^)^ and (B5) • Recalling that Cp fw -C 3 and C 2 ^ 
inserting these expressions for strain components in equation (C4), and 
integrating over the area of the plate, the extenslonal strain energy 
becomes; 

Vep = + «3est^o^ + ^^st^) (C5) 




•i^re with the ■boundary condition by' =0 at y = 0 and y = b 


Mg = + Oi^ + Oa^ + 03 ^ + - V®) -t- SO 5 ® + sa^ + aos^ + 2Bq^ + + Hg^ + H 3 ® + + hb^ + 

+ 0^2 + <^323 + aOjEi^ + 2GgHg + SO^Hy + 20gH^ + 'j[^0*h. SQj^O^ + HqH]^ + 2Hi^ +_ 

v^GqHi + O^Hq + 2 Gj^^ + + -^[ 0(503 - G^Oa ‘ ^<^^6 ^ ‘ ^ 5®6 ■*■ * 

v(aoH2 + GgHg - G^Eg - GgH^ + 20i^a^ + 20^ - EG^g - 80^5^ + ^joaG3 + 20g0^ + HgH^ + SE^ + 
v( 02 H^ ■*■ O 5 H 0 + £G^j + EOjHg^ - -;^|a 3 ag + HyEg + V ^O^g + " ]^[p5*^ "'' ^5^ '*' ''(’^8 ''■ ^^^ 8 ^)] ■*■ 

^Ifl^ + Jg2 + J 32 + ^ J2(5Ji - 3 J 3 )] + + v)® - CiCg(l - 2v - 3 v 2) + C22(l - V + 2v2)| + 

|;[ci(l- V) J.[^^(i ■ v) - ^Ho ^ C2) (4 ^ pg) -ac^Gi - 

^ ^ 2 ) (V ^ ^ Ga) (1 . - 02(X - P^^Ga - 

■ + p^) - Cs(l - P?|h 2 -<-- ^f^ - ~^jg |(ci + 9p2) - 8C^G3 - 

g*(l + p2j Li -i + 9p2j Li -I 


4 ^ h - ^ - «?l-3 - 




g^d - v) L. , i sv 2^' \ f,~lT gp'(i - v) [Z , / gy iLji 

( 1 1 flh y * 3 ■ 1 • sT!l " ■ .(k • 9s")r 


2ic 


4 + p‘ 

9PV J 


M3 - -(1 - v^[i(3 + f?) + -^] 

Hj^-|(l-v") 




(C6) 


(£■ 
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and vlth tlie 'boimdary condition v = 0 at y = 0 and y « b. 


M2 = i[soo" + Ox® + 02^ + 03^ + JtOlf® + aoj® + 206® + 2O7® + 208^ + ^[aOoHo + OlHi + OgHe + 03H3 + 1 kJ|,EH + 20^5 + 206H6 + 2O7H7 + 2Q0H^ 

'Tl.iC- 

^ "^“3 “cdi “^“5 '' "^“5 "5“4;j ■" ' ""i“2 ^ “0^3 ' “i“2 ■" ‘ ™5“6 

3 + O3S0 - Oj^Hg - OgHj^ + 20i^^ + 20^1j - aOjHg - 20^5)] + ^[^3 + + HgHg + + v(0gH3 + + 20^^ + - 




— fn n. 4. n _n. 4. »i / n_w_ 4. a_n ^ 2 _fn_ 


n- 4. n_n- 4. 


M (n n* 4. 


7*[^ro *T*o 'V'TTJ ■ nr7^j i5,[^“o ’ “5“o ' 'V3“o ' '■tr?^ ■ 16 1:1 


Tl 4. ^ .7 V [t ^ 4. T ^ 4. -T ® X T ^<^^ , 


15« 


■2\^i ■ 


I ^ - 2V - 3V®) + Ce®(l - V + 2v®)] + I (5^6® + P®06(vO^ + Hi,) + |j |ox(l - v) + W . c^Qq - 

^ .,\ _ a(l - v) , P®(1 - v) X ^ 3v li, X a^\ _ _ ^®(l - v)_ f7n_ _ 1 - V f,, x =2'l _ nn"l o X 

2 « Ql'- - V - 1 + V - “2)“0 ' -^JV ■ - / ^ V— - 

^ |l + ^ - ^ 2(1 - ~ ^^^2 [ex - Cfi)(l + p®) - CgU - P®^H2 + 

- “!l“= - I' - '»)(*•’'=) - *3-3 - - hf)=! 

• (3^ - 7f7)“3'= - »E- (^ ■ 

Hj - - Ajs + + A®(l + p^v)] - -^(l + V,®)) 


> (07) 


Hi, - I 


y 
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The extenslonal strain energy of a stiffener is 


'es 


1 2 
“ 2 ®st^st / f^x ' ) 

'-'0 ^ / y=Tl^ 


(C8) 


at y = Tjh; write 


’ = qp^t^ cos 2cpx + N2) - Cg^j 


(C9) 


where 


% = 2 ^Ci + C2)(cosh 2 (prih - sinh 2cpT]h) - itC2cprih(cosh 2cpTih - sinh 2 qpiih) + 
Gq + &! sin PqpTih + G2 cos 23 qprih + G3 sin 3 Pcpilh (CIO) 


and 


N2 - Gi^ + G^ sin PqjTih + Gg cos 23qpnh + Gy sin 339nt> + G3 cos ^pcpt^h 


(cn) 


Integration of equation (C8) therefore gives 


^es “ ^ ^^"^^st^st^ 


Wo^-|- + N2^] - 2Wg%r 


_ 2 

st ^st 


(C 12 ) 


The tending strain energy of a stiffener is 


Vhs 2 ®stlst ^ 


(CI 3 ) 


since 


2 

(^\x) _ yj ^ ~ cos 2^^n\b + A sin Pcprjh) sin qjx 

integration of equation (C 13 ) gives 


(Cll^) 
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h 6 

where 


= (1 - cos 23cpq‘b + A sin pcpi^h)^ 
The twisting strain energy of a stiffener is 


(CI5) 



(CI6) 


Since 


(w ) = q)^-tWQp|i sin 2P{pT)h + A cos ptpn'^ cos CpK 

'' ^'y=r\h 


integration gives 


Vts = i <P^t2caWo2 ^ 

where 


(CIT) 


N). = 


3^(2 


sin 2P(prjh + A cos Pcprjh)' 


(C18) 
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CASE 2 


CASE I 



("b) Case 1 , plate deflecting around stiffeners; case 2 , system displacing 

as con5)lete unit. 

Figure 1 .- Sketch of reinforced plate. 
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Figure 4 .- Variation of Kmin with li. 
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FIXED END 


a.n torsional rigidity factor 


SIMPLE 

SuS»dRT 


Figure Variation of ^ correeponding to vltii n 

B«Cs denotes "boundary condition," 
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Figure 9«- Con^iarison of transverse wave form, x = X./2. 
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p = BENDING RIGIDITY RATIO 

Figure 11 .- Variation of p corresponding to K^nin vl'tli P. 
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